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Abstract 

We extend a recent global analysis of nucleon parton distributions carried out at next- 
to-leading order (NLO) in the MS scheme to provide distributions in the so-called DIS 
scheme. We pay particular attention to the translation of the heavy quark distributions 
in going from the MS to the DIS scheme. We repeat the global analysis at leading order 
(LO) and comment on the major effects produced when going from LO to NLO. Finally 
we include in the global analysis a freely parameterised form of possible higher twist 
contributions to make an exploratory study of the size of these effects as a function of x. 


1. Introduction 

We recently published a global analysis (MRST)Q [Q] of data in which we extracted quark 
and gluon distributions for the proton. The aim of the analysis was to constrain the partons 
by the present data — on deep inelastic scattering (DIS) at HERA HI |^, from DIS with 
hxed target experiments |§|, ||, |^, |], and from prompt photon production (PPP) |^. In 

addition the Dreh-Yan measurements [^, together with the pp/pn asymmetry in Dreh-Yan 
HT^ [T^ and the asymmetry of the rapidity distributions of the charged lepton from —>■ 


decays at the Tevatron pp collider |]T^ were used in the analysis. Particular attention was given 
to the PPP data and the influence they have on the determination of the gluon at large x. 


To obtain a satisfactory description of the higher energy data some intrinsic transverse 
momentum of the partons needs to be included and by varying the value of (/ct) within an 
acceptable range we arrived at a range of estimates for the gluon distribution at large x. The 
parton distributions corresponding to the extremes of this range were labelled MRST( 5 f |) and 
MRST( 5 f 1) and these together with the set using the ‘central’ gluon, MRST, all used a value 
of as{Mz) = 0.1175. In order to reflect the range of uncertainty in as^M^) we also provided 
sets of partons (with the central gluon choice) where as{M‘^) was varied by ±0.005 and these 
sets were labelled MRST(a 5 t|) and MRST(a 5 J,|). 

Our analysis included next-to-leading order (NLO) corrections which were computed in the 
MS factorisation scheme. A special feature of the analysis was the new treatment of heavy 


flavour production in DIS, using the procedure developed by two of us |^, |^, which describes 
the threshold behaviour correctly and which is consistent with the MS scheme. The resulting 
hve sets of parton distributions were therefore appropriate only to processes evaluated in the 
MS scheme. 

In this paper we extend our previous analysis to confront a variety of issues. In particular 
we present^ (i) parton distributions also at NLO but computed in the so-called DIS scheme ||I8[| , 
(ii) parton distributions resulting from simply a LO analysis and (iii) parton distributions re¬ 
sulting from a NLO analysis which in addition allows for an empirical higher twist contribution. 


2. Partons in the DIS scheme 

The DIS scheme is simply a device for absorbing the one-loop MS coefficients into a re- 
dehnition of the parton distributions so as to exactly preserve the NLO value of F 2 but using 
an apparent LO expression. No new global htting procedure is involved. 


^In this paper we use ‘MRST’ to denote our previous parton analysis Q in the MS factorisation scheme. 

^The FORTRAN code for all the parton sets described in this paper together with those of MRST can be 
obtained from jittp : //durpdg. dur. ac .uk/HEPDATA/PDF , or by contacting W.J.Stirling@durham.ac.uk. In 
addition, because of the slightly complicated expressions involving heavy flavours, the packages for computing 
the structure functions from each set are provided there. A FORTRAN routine to compute the charged current 
structure functions ±2 and xF^ from the MRST parton distributions has also recently been included. 
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The relation between a quark density in the DIS and MS schemes is, at NLO, 

+ (^)cy(.-)0 9^(i/j,Q^) (1) 

where and are the normal MS massless coefficient functions, e.g. 

= d“’(‘) los((l - ■■)!■■) + 8^(1 - - 1- 

The corresponding relation for the gluon is a matter of convention and normally one hxes it to 
maintain the conservation of momentum. 

- 2nf(^')cf‘>(z)»g^(x/z,Q% ( 2 ) 

where T,{x, Q"^) = J2^£i[qi{x, Q^) +qi{x, Q^)], and where Uf is the number of light flavours. The 
attraction of the DIS scheme, of course, is that if to LO the F 2 structure function is given by 

F 2 {x, Q"^) =J2x [a 2 ,i qi{x, + 02 ,* qi{x, Q^) ] (3) 

i 

then the same expression holds at NLO provided qi = qf^^. However in this scheme it is only 
the F 2 structure function that has no explicit higher order perturbative corrections. Thus, if 
at LO xF^ is given by 


xFsix, = XI ®3,i qi{x, Q‘^) + a^,i qi{x, ], 


( 4 ) 


then at NLO we have 
xF 3 {x,Q^)= 5(1 - ^(1+ 2 :) 


q?^^{x/z,Q^) +a3,i q^^^{x/z,Q^) ]. (5) 


This is the prescription for light favours. For the case of heavy quarks then in some scheme 
which uses MS evolution for the parton distributions, such as that in MRST, the coefficient 
functions for the heavy quarks are mass dependent. Hence we have to decide how best to dehne 
a DIS scheme for massive quarks. The approach we have adopted is to make exactly the same 
change of dehnitions for the partons as above, i.e. using massless coefficient functions only. 
Hence the heavy quarks in this DIS scheme evolve exactly like the light quarks in DIS scheme 
(in the MS scheme the evolution starts from zero at /i^ = rnjj, so from Eq. (^ we see that in the 
DIS scheme it starts from a nonzero value dependent on the gluon distribution). However, the 
DIS scheme coefficient functions for heavy quarks are nontrivial, unlike the case for the light 
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quarks. Expressing first the charm structure function in terms of the MS (MRST) charm 
quark we have 
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and then using Eq. (|^ to dehne the DIS charm quark in terms of the MS (MRST) charm quark 
distribution we hud 
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(7) 

where both coefficient functions depend on rric (and where in principle there is a NLO quark 
coefficient function which we omit due to its insignihcance in practice). All the relevant coef- 
hcient functions can be found in [^. The gluon coefficient function is actually discontinuous 
at = ml, countering the fact that the charm evolution starts from a nonzero value, to give 
a continuous structure function. The hnal term above, involving the double convolution of two 
coefficient functions and the gluon distribution, is potentially rather complicated, especially 
since C^l is dehned in terms of a convolution itself. However, in practice we hnd that using 
Cfl(z) = (1 - - z) in this term alone (making the gluon coefficient function 

continuous) gives an extremely good approximation to the true result; the contribution from 
the nonzero charm distribution for immediately above being negligible. In the limit 


d)/ 


Q'^/rril —>• cx), when —>• 5(1 — ^) and C 2 ^g{z) —>• ^ 2 ^ 3 ( 2 ^), the usual trivial relationship 

between the structure function and parton distribution in the DIS scheme is regained. 

We could alternatively have dehned the heavy quark DIS scheme by demanding this trivial 
relationship for coefficient functions. This would lead to very complicated, mass-dependent 
splitting functions in the DIS scheme. Hence we prefer to keep all the mass dependence in the 
coefficient functions, as in the MS scheme, and let the partons be the usual DIS partons. 

We are therefore able to provide the parton distributions in the DIS scheme analogous to 
those in MRST in the MS scheme [§, labelled MRSTDIS, MRSTDIS(c/ T), MRSTDIS(c/ i), 
MRSTDIS(as TT) and MRSTDIS(as ii)- 


3. Leading order parton distributions 

For some purposes, e.g. Monte-Carlo simulation programs, partons distributions which 
attempt to describe the data at the leading order level are preferable. To obtain such partons 
we repeat the global analysis at LO. This means that the partons evolve only via the LO 
DGLAP equations and each process is expressed in terms of the partons via LO formulas. The 
starting distributions at Qo = 1 GeV^ have the same functional form as in the NLO analysis 

i- 

The heavy flavour contributions to the structure function F 2 are computed to LO which 
means that for we have 

F^{x, Q^) = ^ m^/Q^) 0 g{x/z, Q^), ( 8 ) 
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while for 


F^{x,Q^) = F^{x,Q^ = ml) + -Cf\z,ml/Q‘^) ® c{x/z^Q'^). 


(9) 


Again, the coefficient functions in Eqs. ( 1 , 1 ) can be found in 


Tel. 


There is a signihcant difference in the sizes of the LO and NLO gluon distributions at large x 
which reflect the importance of the NLO corrections to the prompt photon production process 
(PPP) in that region. Thus to get the same cross section the LO gluon has to be typically 
greater than the NLO gluon by about 40-50% for x = 0.3 — 0.45. The LO gluon is also larger 
than the NLO distribution at small x, this being required for a good description of the small-x 
HERA data, and reflecting the importance of the NLO corrections to the quark evolution at 
small X. This leads to a quite different prediction for the longitudinal structure function (which 
has to be taken into account when obtaining the values of F 2 {x,Q‘^) from measurements of 
the cross section). For example the LO value of Fl is twice the NLO value at x = 10“^ and 

= 4 GeV^. This again illustrates the importance of a precise measurement of Fl at small 
X as a test of higher order perturbative corrections. We note also that the Drell-Yan LO cross 
section requires a phenomenological iP-factor of the order of -|-30% in order to get acceptable 
agreement with experiment. 

To obtain a reasonable description of the DIS data we require a larger value of as^M^) 
than in the NLO case, but too large a value spoils the simultaneous description of T 2 , F 2 and 
the PPP data. Using the simple scale choice of jj? = we hnd a value of Alo( 4 flavours) 
= 174 MeV; a satisfactory compromise which implies a value as(Ff§) = 0.125. We dehne our 
‘central’ LO solution - MRSTLO - with this value of and with a gluon constrained by the 
PPP data as before.^ The MRSTLO(q!s' |t) and MRSTL 0 (a 5 jj.) solutions again correspond 
to varying this central value of as(Ff§) by ±0.005. 


Overall, the quality of the LO and NLO hts are comparable, see Table 1, with two excep¬ 
tions. The SLAG data (which cover a region^ of large x and relatively low Q^) strongly prefer 
the NLO corrections — a fact established a long time ago []^. The effect of these corrections 
at large x is equivalent to using a LO description where the value of 05 increases with n, 
the moment of the structure function 0,0 E3- While this effect helps the description of 
the SLAG data, the BGDMS data actually do not favour this trend. Gonsequently, despite 
the relatively large value of as, the absence of NLO corrections can approximately mimic at 
large x a NLO ht with a lower value of as- Thus the BGDMS data are surprisingly well de¬ 
scribed by our LO global ht which can be compared with the MRST(NLO) description in Fig. [^. 


4. Global analysis including higher twist terms 

In MRST we showed that our NLO ht slightly underestimated the slope dF 2 /d\ogQ‘^ for the 
NMG data since the low data tended to he systematically below the ht. Glearly the inclusion 

^The same value (kx) = 0.4 GeV(0.92 GeV) for WA70 (E706-530 GeV) is used as in Ref. |^. 

“^Only data for x < 0.7 are included in these fits - see the discussion in the next section. 
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Data set 

No. of 
data pts 

MRST 

MRSTLO 

HI ep 

221 

164 

159 

ZEUS ep 

204 

269 

288 

BGDMS pp 

174 

248 

171 

NMG pp 

130 

141 

144 

NMG pd 

130 

101 

115 

SLAG ep 

70 

119 

188 

E665 pp 

53 

59 

52 

E665 pd 

53 

61 

62 

GGFR F^^ 

66 

93 

79 

GGFR F^^ 

66 

68 

83 


Table 1: The values for the DIS data included in the (NLO-MS) MRST £t of Ref. [Q] and 
the corresponding LO fit. 


of a negative 1/Q^ contribution is bound to improve the quality of the description in this x 
range 0.02 — 0.1. Higher twist terms have always been expected to play an important role at 
very large x — indeed we found in the past that the STAG data for x > 0.7 were dominated 
by power corrections and for that reason were excluded in MRS leading twist analyses. In those 
analyses, including MRST, we also imposed a lower cut on all data of 10 GeV^ to reduce 
the effect of unknown higher twist contributions. Now that we are allowing such terms, in this 
section we relax the constraint so that only data for < 4 GeV^ are removed. In addition 
we have lowered the cut on data included in the £t from 2 to 1.2 GeV^. 

At very small x we are also interested in examining whether 1/Q^ corrections may be im¬ 
portant. Some attention has been paid to the observation that at HERA the slope dF 2 /d\ogQ‘^ 
appears to ‘flatten’ off around x = 10““^ in contrast to the naive DGLAP expectation. 


In MRST this was attributed to a ‘valence-like’ behaviour of the gluon at the starting scale 
Ql = 1 GeV^. That is for x below 10“^ the gluon is suppressed and consequently, since for 
the HERA data x and are strongly correlated, this leads to a leveling off of the slope at 
very low values of x and Q^. However this ‘valence-like’ behaviour of the gluon may be an 
artifact, reflecting some dynamics other than DGLAP, and one candidate is a positive higher 
twist contribution which is relevant only at small x. 


We assume a very simple parameterisation of the higher twist contribution to the DIS 
structure function, 

F^^{x, Q^) = F^^{x, > ( 10 ) 

where the leading twist NLO structure function F^"^ is treated exactly as in MRST. Apply¬ 
ing the same overall parameterisation independent of target or beam is probably an over- 
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simplification but may have some justification if the higher twist contribution can be derived 
through renormalon dynamics [p^] . 

We parameterise the coefficient D 2 {x) by a constant over different bins in x chosen to 
emphasize aspects of different datasets. We perform hts using the ‘central’ gluon type of 
solution. In general we hnd that for x < 0.5 the resulting correction is small and negative but 
beyond 0.6 large and positive. The values of D 2 {x) obtained from the £t in each x bin are 
shown in Table 2. 


X 

D2{x) (GeV^) 

0 - 0.0005 

0.0147 

0.0005 - 0.005 

0.0217 

0.005 - 0.01 

-0.0299 

0.01-0.06 

-0.0382 

0.06 - 0.1 

-0.0335 

0.1 -0.2 

-0.121 

0.2 - 0.3 

-0.190 

0.3 - 0.4 

-0.242 

0.4 - 0.5 

-0.141 

0.5 - 0.6 

0.248 

0.6 - 0.7 

1.458 

O 

1 

o 

bo 

4.838 

0.8 -0.9 

16.06 


Table 2: The values of the higher twist coefficient D 2 {x) of Eq. (^) (in GeV^) versus x. 


Looking at the overall improvement of the resulting fit we note no real difference in the 
comparison with the HERA data, indeed the higher twist contributions chosen for x < 0.01 
are very small. We have tried hts where the starting gluon at Qg = 1 GeV^ was forced to 
be hat, or even singular, as a: —0 to see if the ehect of the valence-like gluon could be 
described instead by a significant positive higher twist term. In each case the quality of the 
resulting ht was (far) worse. So there appears to be no preference for a description in terms of a 
‘conventional’ gluon at low with additional power corrections, at least within our admittedly 
simple parameterisation. 

At intermediate x values (0.01 < x < 0.5) there is dehnitely a preference for some negative 
1/Q^ contribution as we expect from the NMG data. With the relaxing of the cut, the 
number of NMG Tf points increases from 130 to 155 but the stays close to 140 with the 
inclusion of the higher twist term. The most signihcant improvement in the quality of the ht 
is for the data for x > 0.1 where the higher twist corrections are largest. The most dramatic 
diherence is the description of the SLAG data. In Fig. ^we illustrate the improvement at very 
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large x due to the addition of the arbitrary higher twist contribution. From this hgure and from 
the values of the coefficient D 2 {x) in Table 2 we see that the very large x SLAG data require a 
large positive 1/Q^ correction. We have not attempted to separate the part arising from target 
mass corrections but a recent analysis suggests a significant fraction can be accounted for 
in this way. 

We may regard the difference between this new set of partons - MRST(HT) - and MRST 
as one measure of uncertainty in our parton sets. In Fig. we show the ratio of the u and d 
partons at = 10 GeV^ between the two sets — the differences being much smaller for other 
partons. We see that the differences remain less than about 1-2% except where the individual 
distributions start to become really quite small. At large x, x > 0.6, the large differences 
for the valence distributions are not surprising in view of the large higher twist correction as 
a: —> 1. As increases the ratios remain rather constant so that Fig. ^ is a reliable indication of 
the uncertainty on the u and d partons at all arising from possible higher twist contributions. 


5. Conclusions 


We have presented an extension of our previous global analysis to provide alternative parton 
distributions. In particular we have obtained the analogues of the partons of MRST but (a) in 
the DIS NLO scheme as opposed to the MS scheme and (b) by repeating the global analysis at 
LO. In this latter case the partons are consistent with all processes considered being evaluated 
to LO. Finally, we present a set of partons obtained from a global analysis in which an empirical 
universal higher twist contribution is included, which is freely parameterised as a function of 
X. This empirical higher twist component is found to be surprisingly small in the HERA small 
X domain. On the other hand it is interesting to see that the higher twist fit is similar at high 
X to that expected from a renormalon approach for the nonsinglet structure function |2^ . 
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Figure 2: Comparison of very large x data from BCDMS and SLAG with the standard NLO 
fit of MRST and the higher twist £t MRST(HT). 







